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Abstract 

We develop the theory of conformal blocks in CFT^ expressing them as power series 
with Gegenbauer polynomial coefficients. Such series have a clear physical meaning 
when the conformal block is analyzed in radial quantization: individual terms describe 
contributions of descendants of a given spin. Convergence of these series can be 
optimized by a judicious choice of the radial quantization origin. We argue that the 
best choice is to insert the operators symmetrically. We analyze in detail the resulting 
"/9-series" and show that it converges much more rapidly than for the commonly used 
variable z. Wc discuss how these conformal block representations can be used in the 
conformal bootstrap. In particular, we use them to derive analytically some bootstrap 
bounds whose existence was previously found numerically. 
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1 Introduction 

Recent years have seen a revival of the bootstrap approach to Conformal Field Theory 
(CFT) in higher dimensions [1-14]. Recall that conformal bootstrap aims to control a CFT 
by imposing the associativity constraint on the operator algebra. In practice this is done by 
taking a correlation function of four primary operators^, expanding it into conformal partial 
waves, and demanding that the different channels agree (see Fig. 1). A central role in this 
program is played by conformal blocks — functions of the cross ratios obtained by stripping 
the conformal partial waves from the trivial x-dependent factors. 

The theory of conformal blocks in dimension d ^ 3 was started in the 70's [15-18], with 

^These are called quasi-primaries in d — 2 dimensional CFT. 
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many recent valuable contributions [19-22, 10, 23, 24].^ Especially the explicit expressions 
for even d found by Dolan and Osborn [19, 20] were instrumental for the first practical 
applications of the bootstrap. In general d, an approach to evaluating the conformal blocks 
and their derivatives was developed last year in [10], and applied in the bootstrap analysis 
of the 3(i Ising model. 



Figure 1: Any CFT is characterized by conformal data — primary operator dimensions and 
spins (Ai,/j) and the OPE coefficients fijk- Using the OPE, the four point functions can be 
expanded into conformal partial waves, fixed by conformal symmetry in terms of the operator 
quantum numbers, times the products of the OPE coefficients. That the different expansions 
agree is a nontrivial constraint on the conformal data. 

In spite of these advances, conformal blocks remain rather mysterious special functions. 
The purpose of this paper is to demystify them via a concrete and economical approach. We 
will base our considerations on the fact that the conformal blocks are, first and foremost, 
sums of contributions of radial quantization states to a matrix element computing a four 
point function. This point of view is standard in the 2d CFT literature [26]. Recently 
[27], it proved useful in general d to study the convergence rate of the conformal block 
decomposition. Here we develop it to its logical end. 

The paper is organized as follows. In section 2 we introduce the radial quantization 
representation of conformal blocks. Here we consider the four point function in the frame 
where two points are fixed at and oo. The expansion parameter in this frame coincides 
with the Dolan-Osborn variable z. The heart of the paper is section 3, where we switch to 
a different coordinate p, which corresponds to the frame with points inserted symmetrically 
with respect to the origin. We demonstrate the advantages of this frame for evaluating the 
conformal blocks: the expansion parameter p is smaller than z\ the p-series converges every- 
where where the block is expected to be regular; its coefficients are bounded independently 
of A and /. The last two properties are not true for the 2;-series expansions. 

We foresee that the p-series representations of conformal blocks will find many applica- 
tions in the conformal bootstrap program; we outline some in section 4. An especially neat 
application is the "toy bootstrap equation" (section 4.2), by means of which it is possible 
at last to get analytic understanding of why the methods of [1] were successful in producing 
upper bounds on the operator spectrum. 

^In d~ 2 dimensional CFT, one distinguishes the "small" conformal blocks defined by summing over the 
SL{2,C) descendants and the "big" blocks defined by summing over the Virasoro descendants. Although 
our focus is on higher dimensions, the blocks considered here reduce to the "small" blocks in d = 2. They 
can also be viewed as the c — >■ oo limits of the "big" blocks [25]. 
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We conclude in section 5. Appendix A contains the proof of boundedness of the p-series 
expansion coefficients. 

2 Conformal blocks in the Dolan-Osborn coordinates 
2.1 General structure 

For simphcity we will focus on the Euclidean-space correlator of four identical scalar pri- 
maries.^ By conformal invariance it has the form 

where g{u,v) is a function of the conformally invariant cross ratios 

— (^12^34)/(^13^24)> — (^14^23)/ (^13^24)- (2-2) 

The partial wave decomposition of this correlator takes the form: 

giu,v) = J2foGo{u,v), (2.3) 
o 

where the Go{u, v) are the conformal blocks of the primary operators appearing in the 0x0 
OPE and fo are their OPE coefficients. The function g{u, v) computed from this expansion 
must satisfy the crossing symmetry equation 

v^^g{u,v) = u^^g{v,u), (2.4) 

which imposes constraints on the dimensions, spins, and OPE coefficients fo of the ex- 
changed operators. However, our main interest here is not in how to extract these constraints 
(this will be briefly discussed in section 4), but in the conformal blocks themselves. 

Starting from the work of Dolan and Osborn [19, 20], it has become customary to express 
conformal blocks by changing coordinates from u,v to z and z = z*: 

u = zz, V = {I- z){l-z) . (2.5) 

The geometrical meaning of the new variables is made clear by assigning three points to 
0, 1, oo as in Fig. 2. The complex z is the usual coordinate used in c? = 2 dimensional CFT, 
but its utility for general d is not a priori obvious. Refs. [19, 20] discovered that conformal 
blocks in = 4 and in all even dimensions take particularly simple expressions in these 
coordinates. Here we will work with any d, even or odd."^ 

■^The generalization to non-identical scalars is straightforward, and we will comment on the non-scalar 
case in section 5. 

■*In fact, as we will see below, conformal blocks depend analytically on d. The conformal bootstrap 
equation with analytically continued blocks can be formally considered for any d. It can be taken as a 
nonperturbative definition of CFT in fractional dimensions [28]. 
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Figure 2: By conformal symmetry, three operators can be put at x\ = 0, X3 = (1, 0, . . . , 0), 
X4 — )• 00, with the fourth point X2 somewhere in the (12)-plane parametrized by the complex 
coordinate z. 



To avoid possible misunderstanding, we should stress that although we parametrize the 
conformal blocks by a complex variable we never use complex analysis. Only in the 2d 
case do the conformal blocks factorize as a holomorphic times antiholomorphic function. 
For general d considered here, we will treat conformal blocks as smooth real functions in 
the z plane; see the end of this section and section 3.2 for more details. 




Figure 3: Using a Weyl transformation, the configuration in Fig. 2 is mapped onto a cylinder 
matrix element with operators inserted as shown. 



To exhibit the general structure of conformal blocks in the z^ z variables, let's use radial 
quantization. It is convenient, although not strictly necessary, to do a Weyl transformation 
which maps the CFT from to the cylinder M x S'^~^ . In polar coordinates, the mapping 
is simply (r, n) — t- (r, n), with the cylinder time r = logr. The flat space four point function 
with points assigned as in Fig. 3 then maps to the cylinder matrix element 

(</)|0(r3,n3)0(r2,n2)|0). (2.6) 

The operators inserted at zero and infinity map to the radial quantization in- and out-states 
|(/)) and ((/)|. The other two insertions are at the cylinder times T2 = log \z\ and rs = 0. We 
keep both unit vectors 112 and explicit for future use, but the only rotationally invariant 
parameter is their scalar product 

n2 • n3 = cos 9, 9 = arg z . (2.7) 
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The next step is to express (2.6) by inserting a complete basis of energy eigenstates on 
S'^~^. This gives^ 

g{u,v) = J2\^f{mO,n,)\E) {E\<j){0,n2M) (2.8) 



E 

where we took into account that propagating a state of energy^ E for the Euchdean time 
distance T3 — T2 will give rise to the factor 

g-i?(r3-r2) _ ^2.9) 

The exchanged states on the sphere are in one-to-one correspondence with the local 
operators appearing in the OPE 0x0. For the moment we do not distinguish between the 
primary and descendant states. Every state will come in a multiplet of SO{d). In fact, only 
symmetric traceless tensor multiplets of spin j ^ can couple for the considered correlator.^ 
The right matrix element 

(E,{/ii,...,/i,}|0(O,n2)|0) (2.10) 

must be a rank-j symmetric traceless tensor constructed out of the vector 112, which is fixed 
up to a constant: 

nfnf...n^^ -traces. (2.11) 

Analogously, the left matrix element is fixed up to a constant, and so a general term in (2.8) 
will be proportional to^ 

(nf nf . . . ni^' - traces)(nf nf . . . n|^^' - traces) oc CJ{n2 ■ ng), u = d/2 - 1, (2.12) 

where Cj are the Gegenbauer polynomials. For the integer dimensions of interest d = 2,3, 4, 
they take the form 

lim z/-^a^(cos 9) = - cos(j 9) (j ^ 1), 

Cy \cos 9) = Pj {cos 9), (2.13) 

sm 9 

In particular, for = 3 we get the Legendre polynomials. 

We conclude that the function g{u,v) appearing in the four point function (2.1) must 
have an expansion of the form: 



g{u,v) 



1 + I]Pi^jl-s|''CJ(cos^) , Ve,j > 0. (2.14) 



^The extra factors in the RHS of this formula which follow from the denominator of Eq. (2.1), from the 
Weyl transformation of operators, and from acting with exp(±i?T) on the in- and out-states when shifting 
the operator insertion times to zero, cancel each other; see [27] for a more detailed derivation. 

^The presence of the Casimir energy on the sphere, nonzero for even d, can be ignored here. This is 
because we are discussing correlation functions and not the partition function, and so the relevant energy 
is the one defined subtracting the energy of the ground state. 

^One cannot construct an antisymmetric tensor out of a single vector n, and so the corresponding matrix 
elements necessarily vanish. 

^This contraction formula follows from the theory of spherical harmonics and harmonic polynomials, see 
[29], Section 11.2, Lemma 1, and [30], Chapter 4. 
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where the sum is over all local operators of dimension E and spin j appearing in the OPE 
(j) X (f). Although the coefficients pE,j are left undetermined by this argument, we do know 
that they must be non-negative. This is because for 112 = the configuration in Fig. 3 
becomes reflection-positive. The matrix elements in (2.8) are then complex conjugates of 
each other. 

The appearance of Gegenbauer polynomials in this result is not surprising, as they 
already arise in the theory of angular momentum in quantum mechanics. When two spinless 
particles scatter through a spin-j resonance, it is well known that the amplitude is given by 
the Legendre polynomial of the scattering angle (see Fig. 4). 




Figure 4: Elastic center-of-mass scattering of two scalar particles. When a spin-j resonance 
dominates the scattering process, the amplitude is proportional to Pj(cos6). 



Consider now a particular primary operator O of dimension A and spin / occurring in 
the (j) X (f) OPE. If we restrict the sum (2.14) to its conformal multiplet, it must represent 
the conformal block of O. The conformal multiplet will have descendants of integer-spaced 
dimensions A + n with spins at level n taking values^ 

j = I + n,l + n — 2, . . . , max(Z — n,l + n mod 2) . (2-15) 

Moreover, the coefficients pE,j within one conformal multiplet are not independent, since 
the matrix elements for the descendants will be all proportional to the basic OPE coefficient 
fo- We conclude that the conformal block must have the following expansion: 

^ ^-^ Q(cos^) 

n=0 j ) 

where the positive coefficients Anj are some universal functions of A, /, and d that are fixed 
by conformal symmetry. We normalize the total conformal block by the condition Aq,; = l.^'' 

The Gegenbauer normalization factors, 

q(l) = (2z.),/j! , (2.17) 

are included in (2.16) for later convenience and also to ensure a smooth u limit for 
d = 2. Here, (x)„ = r(x + n)/r{x) are the Pochhammer symbols. 



'For the short representations some of these spins wiU not be there. 



^"This normahzation relates to the one used by Dolan and Osborn in [19, 20] as G^'^° 



(-2)'(;.),/(2i.),G"h- 
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The formula (2.16) is the main result of this section. It should be noted that Ref. [20] 
already used an expansion of conformal blocks into Gegenbauer polynomials, because they 
turn out to form a convenient basis for solving the Casimir differential equation recursively 
(see the next section). Ref. [13], Eq. (78), observed that in any number of dimensions 
conformal blocks can be expanded in cos(j 6) with positive coefficients. For d = 2 our result 
says the same, although for general d our conclusion is stronger. To obtain their result, 
one runs the above argument classifying states into multiplets with respect to the SO {2) 
subgroup of SO{d) acting in the (12)-plane. In particular, the Gegenbauer polynomials for 
any z/ ^ have positive expansions in cos{j' 6), j' ^ j. 

The region of convergence of the expansion (2.16) will be limited to \z\ < 1, which is the 
condition for the operators 02 and 03 in (2.6) to be time-ordered on the cylinder. However, 
the actual domain X of regularity of the conformal block as a function of z is larger; it is 
given by the complex plane minus the (1, +oo) cut along the real axis: 

X = C\(l,+cx)). (2.18) 

Everywhere in this region the blocks will be real analytic, except at 2; = because of the 
1^1^ factor. For every point in X one can find a sphere which separates Xi and X2 from 
and X4. Choosing the center of this sphere as a radial quantization origin, one can prove 
the regularity of the conformal block for such z. The blocks will be singular on the cut, 
because the separating sphere jumps when z crosses it. In section 3 below we will construct 
expansions convergent in the full region X. But first we would like to study the coefficients 
of the expansion (2.16). 

2.2 Expansion coefficients from the Casimir equation 

We would like to compute the coefficients Anj in (2.16). In principle, this can be done 
following the radial quantization method to its logical end: imposing the constraints of 
conformal invariance in the OPE and evaluating the norms of the descendants. The example 
of scalar exchanged primaries and their first two descendant levels was considered in [27]. 
However, it is far more efficient to use a different method first proposed in [20]. 

The idea is that the conformal block satisfies an eigenvalue equation of the form 

VGa,i{^^,v) = Ca,iGaAu,v), Ca,i = A{A-d) + l{l + d-2), (2.19) 

where P is a second-order partial differential operator. To get it, one acts on the four point 
function with the combination of conformal group generator (in the SO{d + 1, 1) notation) 

l{L^il + L^ll)iLi^^^^ + Li^^^% (2.20) 

where the generators L^'^^ acts on the operator inserted at Xj. By conformal invariance 
of the OPE, this combination can be pushed through to act as the quadratic Casimir on 
the operators appearing in the OPE (j){xi) x 0(0:2) • All terms within a given conformal 
family will have the same Casimir eigenvalue Ca,i- This gives a differential equation for the 
conformal partial wave. 
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In the z^z coordinates the operator D takes the form [20] 

77 

\V = [z\l - z)dl - z^d,] + [z\l - z)dl - z^d,] + 2z/— ^ [(1 - z)d, - (1 - z)d,] . (2.21) 

z z 

For our purposes it will be convenient to express it in the coordinates 

5 = 1^1, ^ = cos9={z + z)/{2\z\). (2.22) 

We find: 

V = Vo + Vi, 

Vo = s^dl + {2u +l)[idi:-s d,] - (1 - e) dl (2.23) 

V, = s [-^s'di +2{i-e)s dsd^ -^sd,- {2u + e) 5? + e (1 - e) ^i] • 

The terms are grouped in such a way that preserves homogeneity in s while T>i increases 
it by 1. 



We now apply this operator to (2.16), which we write as 

{2v 



n=0 j 

Using the properties of Gegenbauer polynomials, it is easy to see that Vej are eigenfunctions 
of "Do- The eigenvalue depends on the dimension and spin in the same way as the Casimir: 

VoVej = Ce,jVej . (2.25) 
The Pi also acts simply in this basis: 

I^iVej = -Iej'Pe+ij+i - -/ejVe+ij-i, 

_ {E + jfU + 2u) _ (E-j-2ufj 

2(i + .) ' 2(j + .) • ^'-'^^ 

Applying these formulas, equation (2.19) can be solved order by order in s. We find that 
the coefficients An^j must satisfy the following recursion relation: 

(CA+n,j - C'a,/) ^nj = 7A+n-l,i-l^n-l,J-l + 1 /\+n-l,j+l^n-l,j+\- (2.27) 

Starting from the initial conditions 

Ao,i = 5,1 (2.28) 
this recursion determines all coefficients An.j- One can check that 

CA+n,i - CA,i > (2.29) 

if A satisfies unitarity bounds and j is in the range (2.15). So the coefficients generated by 
the recursion are manifestly positive, in agreement with the previous section. 
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For illustration, here is what the solution at the first two levels looks like: 



_(A + /)(/ + 2i/) (A-/-2I/)/ 
A-id+i, Al l I — 



A2,l+2, A21, A2I- 



A{l + iy) ' A{l + u) ' 
(A + 0(A + / + 2)^{l + 2iy){l + 2iy+l) 

32(A + / + l)(/ + i/)(/ + z/ + l) 
(A + /)(A - / - 2z/)[(A - + 2v) + (A - 2v){v - 1)] 
16(A-z/)(/ + z/ + l)(/ + z/-l) ' 
[A -I- 2u){A -l-2u + 2fl{l - 1) 
32(A -l-2v + + u){l + u -I) ■ 



(2.30) 



Notice that low spins do not require a separate treatment: the coefficients which "do not 
exist", like for / = and A2^i-2 for / = 0, 1 come out automatically zero. This follows 

from the fact that T^Jo = and so (2.26) makes sense also for j = 0. 

The recursion (2.27) has been found previously by Dolan and Osborn [20], Eq. (3.12), 
who arrived at the ansatz (2.24) as the way to diagonalize the homogeneous part of V. 
They were expanding in Jack polynomials symmetric functions in two variables z, z, which 
are identical to our Vej- They also give a closed-form solution of this recursion, Eq. (3.19), 
which is however rather complicated (it involves 4-F3). In practice, it may be faster to 
evaluate the coefficients directly from the recursion. 



2.3 Decoupling of descendants for the leading twist^^ 

One interesting special case where the recursion can be solved easily is for the "leading 
twist" operators O of dimension 

A = Z + c;-2, / = 0,1,2... (2.31) 

In this case we find that at each level, only the maximal allowed spin j = l + n has a nonzero 
coefficient. At the ffist two levels, this can be seen happening in Eq. (2.30). For general n, 
this single nonzero coefficients takes the form:^^ 

{l + uUl + 2u)^ 
n][2l + 21^ )n 



The massive decoupling of descendants implied by this result can be understood as 
follows. The descendants at level n are obtained by acting with n derivatives 

d,,d,,...d,„0. (2.33) 

If a Hi is contracted with an index of O, such a state simply vanishes, because for / ^ 1 the 
dimension (2.31) is the minimal value allowed by the unitarity bound and corresponds to a 

^^This section is independent of the main line of reasoning and can be skipped on the first reading. 
^^For d = 3, this result is agreement with the integral representation in [21], Eq. (6.20). 
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conserved current. If some of the /Xj are contracted with each other, we get a state involving 
d"^ which has spin strictly less than / + n. We should show that such states decouple. Since 
they do not have zero norm, this can only happen via vanishing of the matrix elements in 
(2.8). Equivalently, this means that the following limit of the three point function should 
vanish: 

lim \xi\^^^{cj){xi)cP{x2)dlO{y)) = (2.34) 

Since the three point function {(p4>0) is known explicitly (see e.g. [19]), this is easy to check. 
Sending xi — )■ cxd, 0:2 — )■ 0, the three point function becomes 

{<P\mO,,...M) = My,. ■ ■■yj\y\''~^^^' - traces) 

°^^Mr--^A«.^pz^- (2-35) 

That the second line takes care of the trace subtractions in the first line (up to a constant 
factor) is obvious: it gives a tensor which has the right scaling in y and is also automatically 
traceless (as well as conserved), due to the fact that the function l/\y\'^~'^ is harmonic in d 
dimensions. For the same reason, this formula implies that 9^-descendants decouple. 

We should stress that the decoupling of ^^-descendants at leading twist is peculiar to 
the kinematic configuration of Fig. 2. In particular, it will not happen when the points are 
inserted symmetrically with respect to the origin, as in the next section. This is because 
Eq. (2.34) is only true in the infinite xi-limit. 



3 Conformal blocks in the p coordinates 

We now wish to analyze the four point function (2.1) in a different, more symmetric, config- 
uration of operator insertions, shown in Fig. 5. Applying a conformal transformation, the 
configuration of Fig. 2 can be mapped to the new one. There is a one-to-one correspondence 
between the complex parameters z and p, fixed by demanding that the cross ratios should 
agree. We find 

z 4p 

P= , <^ Z=- ^-TTT • (3.1) 



The p coordinate was considered in [27], where it was used to give an optimal estimate 
for the convergence rate of the decomposition of a four point function as a sum of conformal 
blocks, Eq. (2.3). Here we will use p to analyze the blocks themselves. As discussed at 
the end of section 2.1, the blocks are expected to be regular in the region X = C\(l, -|-oo). 
The function p{z) maps this region onto the unit disk (see Fig. 6). This suggests that this 
coordinate should be particularly suitable to analyze the blocks. To begin with, conformal 
block representations as power series in p will converge for |p| < 1, which is the full region 
of interest. Other advantages will be discussed below. 

Fig. 7 shows what the configuration of Fig. 5 looks like after the Weyl transformation 
to the cylinder. This picture is similar to Fig. 3 in that both the initial and final state are 
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Figure 5: This more symmetric configuration of operation insertions can be obtained from 
the one in Fig. 2 by a global conformal transformation. 




Figure 6: The p coordinate maps the regularity domain X onto the unit disk. 



characterized by just one unit vector. For this reason the whole discussion of section 2.1 
expressing the exchanges of spin j states in terms of Gegenbauer polynomials goes through 
unchanged. We can therefore state the following analogue of Eqs. (2. 16), (2. 24): conformal 
block of a dimension A, spin / primary will have an expansion: 

oo 

Ga,« = Y.11 ^n,i^A+n,i(r, Tl), 5,,, ^ 0, (3.2) 

n.=0 j 

where 

r = |p|, r7 = cosargp. (3.3) 

The non-negative coefficients i?„ j in this new expansion will of course be different from 
Anj. The spins j at level n will still be subject to the constraint (2.15). However, notice 
that only even spin states can be exchanged since the initial state is symmetric with respect 
to p — 7- ~p}^ We conclude that only even levels n will have nonzero Bnj-^'^ This is unlike 
in (2.24) where all levels have Anj 7^ 0. 

We now turn to the problem of determining the coefficients -Bnj- The first method is 
to convert from the old expansion (2.24) whose coefficients Anj we already know how to 

^■^In fact the exchange 1^2 corresponds to z z/(z — 1), which is equivalent to p ^ —p. 

^^In an analogous expansion for a four point function of non-identical primaries, states of all levels will 
be exchanged. However, if Ai — A2 and A3 = A4, then again only even levels will appear. This is even 
though the exchanged primary may have both even and odd spin in this case. 
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Figure 7: The analogue of Fig. 3 for the new configuration. 



compute. From (3.1), the relevant variables are related by: 



1 + 2r?7 + ' 1 + 2r?7 + 

Substituting into (2.24) and expanding the denominators, we will get a power series of the 
form 



oo 



(3-5) 

n=0 

with Qniv) certain polynomials in rj. To extract -Bn,j, we have to reexpand Qn{v) i^ito the 
basis of Gegenbauers. This will give B^j at level n as a linear combination of A„/ _,/ for 
n' ^ n. 

The second method is to set up an independent recursive procedure for Bnj based on 
the Casimir equation. The operator V in r,7] coordinates takes the form: 

V = Vo + V, (3.6) 

where the homogeneity-preserving part Vq is the same as in (2.23) with s — )■ r, ^ — )■ r^. The 
homogeneity-increasing part is given by 



V = 



1-2^' + ^' ^^rd+ (37) 



1 + r4 - 2r2(2r/2 _ 1) I - 
Its action in the Vej basis will look like 

^^^.^■ = - E ErS'^'^i^w- (3.8) 

n=2,4,... / 

The series is over positive even n, since only such powers of r occur in the expansion of V. 

The dependence of the F coefficients on j' is found with the help of the following identities 
involving the Gegenbauer polynomials (the radial dependence of VE,j is not important here): 

i2r]^ - l)VE,j = a-VEj-2 + a^jVEj + a+P^j+a , 
2r]il - 7f)d^VE,j = b-VE,j-2 + h]VE,o + &^'^i?,i+2 , (3-9) 
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where 



a,- = TTT-. rr-. 7T, dn = T-. TTT-. 7T , CtJ ~ 



' 2(j + z/)(j + z/-l)' ^ + + + ^ 2(j + z/ + l)(j + z/) ' 

^ j(j-l)(j + 2^) ^ J(j + 2z/)z/ ^+ ^ (j + 2z/ + l)(j + 2z/)j 

^' 2(j + z/)(j + z/-l)' + + + i 2(j + z/+l)(j + z/) ■ 

(3.10) 

For example, for n = 2 we get 

The recursion relation for the takes the form: 

(Ca+uj - Ca,i) Bnj = ^ y^r^;];^;-^.,i?„/j/ , (3.12) 

n'=0,2,...n-2 j' 

At level we have the initial condition: 

Boj = k5ji. (3.13) 
We will set = 1, keeping in mind that the normalization of section 2 would correspond to 

To find the Bnj up to level A^, one needs first to compute the coefficients F^^"'-' for 
n ^ N. For example, Eq. (3.11) is sufficient to find the solution for level 2: 

Z(/ - 1)(A -/ - 2z/) Az/(z/ - 1) + (A - l)/(/ + 2z/) 

^2,1-2 = W7i~. 7vn~. — ^7^^^ ; — — r~rT5 = v- 



2(/ + z/- l)(/ + z/)(A-/-2z/ + l)' (A-z/)(/ + z/ + l)(/ + z/-l)' 

R (A + /)(/ + 2z/)(/ + 2z/ + l) 

- 2(A + / + l)(/ + ^)(/ + z. + l)- ^^-^^^ 



3.1 Comparison between the z and /? expansions 

We have presented two ways to expand the conformal blocks: the "z-series" (2.24) and the 
"p-series" (3.2). We will now argue that the second expansion is more efficient, in the sense 
that it converges more rapidly and fewer terms need to be evaluated in order to get a good 
approximation. This happens because of the better choice of the expansion parameter and 
the better asymptotic behavior of the series coefficients. 

Let us start with the expansion parameters. The interesting range for the p coordinate is 
the unit disk |p| < 1. The p-series will converge absolutely, everywhere in this disk. For any 
e > the convergence will be uniform for |p| < 1 — e. To prove this statement, consider first 
p = 1 — e real. For such p all terms in the series are positive, and so the series must converge: 
a divergence here would mean a physical singularity for the conformal block, while as we 
discussed all such singularities are confined to |p| = 1. Convergence everywhere else in the 
disk will be only better. Physically, this follows by the Cauchy inequality: every term in the 
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series is the product of two matrix elements, which for real p become Hermitian conjugates 
of each other. Formally, this is because the Gegenbauer polynomials with our normalization 
are less than one in the absolute value on the interval [—1, 1]. 

The same argument can be used to show that the z-series will converge absolutely in 
the disk \z\ < 1. As we discussed, this does not even cover the full regularity region of the 
conformal blocks. Moreover, from the second Eq. (3.1) we have 

\z{p)/p\>l (|p|<l). (3.15) 

So even in the region where both series converge, the p-series will always have a strictly 
smaller expansion parameter. 

An additional bonus appears when considering conformal blocks for equal external 
dimensions. As we have seen, in this case the p-series involves only even levels. So, the 
effective expansion parameter becomes p^. In conformal bootstrap applications, one usually 
uses conformal blocks evaluated near z = 1/2, which would correspond to p = 3— 2v^ ^ 0.17 
and p2 ^ 0.03. 

Let us now examine the expansion coefficients. We are interested in their asymptotic 
behavior when A or / become large. In the large A limit the coefficients Anj at level n grow 
as 

= 0(A"). (3.16) 

For n = 1,2 this can be seen in Eqs. (2.30). The reason for this growth is that the operator 
Vi is second order in dg. Because of this the coefficients in (2.26) are 0{E'^). On the 
other hand the factor in the RHS of the recursion relation: 

CA+n,j - Ca,/ = 2nA + n{n-d)+3{3 + d-2)~l{l + d-2) (3.17) 

increases only linearly in A. So, going up one level in n, the coefficients A^j gain one power 
in A. 

Turning to the second expansion, we encounter a crucial difference. Unlike "Di, the 
operator T) in (3.6) is only first order in dr- So the coefficients F entering the Bn^j recursion 
grow only linearly in E, and this growth cancels when dividing by (3.17). Contrary to the 
previous case, going one level up in the recursion relation does not increase the leading 
power of A. We conclude that the coefficients Bnj remain bounded in the large A limit. 
For n = 2 this is illustrated by Eqs. (3.14). 

Keeping more careful track of the size of the relevant factors, one can show the following 
sharper statement (see appendix A). Each coefficient Bnj is uniformly bounded in the full 
range of A and / allowed by the unitarity bounds, with the bound depending only on the 
level n and on d: 

max^maxS„ jj ^b{n,d). (3.18) 

The region close to the free scalar limit / = 0, A — )■ i/ is understood excluded when taking 
the maximum. As is well known, the scalar conformal block becomes singular in this limit. 
Physically this is due to the fact that the free scalar must be decoupled from everything 
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else. In our representation, the singularity first shows up in the coefficient -82,0 ~ (A — z/) ^, 
see Eq. (3.14), and then feeds into higher levels. 

Let us discuss a bit how the coefficients Bnj grow with n. This growth is related with 
the the behavior of the conformal block for real p — )■ 1. It can be shown using the results of 
Ref. [10]^^ that in this limit the conformal block has a power-like singularity of the form: 

GA,z(real p ^ 1) ~ (^log ^ for d = 2^ . (3.19) 

On the other hand, the p-series representation for real p > takes the form: 

00 

G'A,/(real p > 0) = /3„p^+^ Pn = Yl ^n,, • (3.20) 

n=0 j 

From compatibility with the p — )■ 1 asymptotics, we conclude that for any A and / the sum 
of the coefficients at level n behaves asymptotically as 

(3n ~ n'^'^ (n^oo). (3.21) 

It would be interesting to know how the ratio /3„/n'^~'^ behaves for small and intermediate n. 
The simplest possibility which accommodates both (3.18) and (3.21) is that /3„ ^ c{d)n'^~^ 
for all n, A and /. However, further study is needed to check this hypothesis. 

To finish this section, we would like to demonstrate how the highlighted differences 
between the z- and p-series can be seen in the explicit expressions for the conformal blocks 
available for even d. These expressions [17, 19-21] are written in terms of the functions 

ka{z) = z''/\Fi{a/2, a/2; a; z) (3.22) 

with a = A + / and A — / — 2z/. For large A, the coefficient in the expansion of the 2-F1 
grows as a" ~ A". This is the same growth as in (3.16). However, when the p variable is 
used, the function ka can be transformed using a hypergeometric identity 

A;j4p/(l + p)2] = (4p)'^/Vi(l/2,a/2;(a + l)/2;p2). (3.23) 
As advertised, this is a function of p^ and the expansion coefficients do not grow with a. 



3.2 Relation to Zamolodchikov's uniformizing variable 

We would like to briefly mention a similarity between the change of variables from z to p 
advocated here, and the one proposed long ago by Al. Zamolodchikov [31] in the study of 
2d "big" conformal blocks. His variable is given by 

q = e'^\ T = iK{l - z)/K{z) , (3.24) 

^^This follows for ^ = from the explicit 3F2 representation on the real line, Eq. (4.10) of [10], and 
remains valid for Z ^ 1 by the recursions in Appendix A of [10]. 

^^This section is independent of the main line of reasoning and can be skipped on the first reading. 
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where K{z) is a complete elliptic integral of the first kind 




2 J, 



(3.25) 



The variable r takes values of the upper half plane, parametrizing the universal covering 
of the Riemann sphere with three punctures 0,1, and oo. The point is that in the 2d case, 
conformal blocks factorize as J^{z)J-'{z) where J^{z) is holomorphic in the complex plane 
with branch points at the punctures. So it is natural to view it as an analytic function on 
the universal covering space. The conformal blocks are then given as power series in q. Since 
|g| < 1 in the upper half plane, these series converge everywhere where J^{z) is analytic, 
while power series representations in z converge only for \z\ < 1. 

In the 2d case, the variable g is a more efficient expansion parameter than our variable 
p. For example, the conformal block regularity domain X = C\(l,+oo) is mapped on a 
subset of the complex plane located strictly inside the unit disk, see Fig. 8. For general d, 
when holomorphic structure is absent, the variable p is probably best possible. 



Figure 8: The regularity domain X is mapped in the q plane onto this UFO-shaped region 
inside the unit disk. For example, the point z = 1/2 is mapped to g ~ 0.043. 



4 Potential applications to the conformal bootstrap 

In the previous section, we introduced a new way to represent the conformal blocks, by 
expanding them in the polar coordinates associated with the complex variable p. Our 
interests in the blocks stems from the role they play in the conformal bootstrap program. 
We believe that our new representation will turn out quite useful in this context. Here we 
will list several ideas, leaving their complete development for the future. 

Most existing applications of the conformal bootstrap to CFT in d ^ 3 dimensions 
followed the following scheme proposed in [1]. Due to the complexities of dealing with the 
full bootstrap system, one focuses on just one equation out of the infinitely many pictured 
in Fig. 1: the one in which all four external states are one and the same, scalar primary 
(j)}'^ This equation is obtained by substituting the conformal block expansion (2.3) into the 

^^The case when the external states are different components of a scalar global symmetry multiplet has 
also been considered [5, 7-9]. 



0.6 r 
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crossing symmetry constraint (2.4) and takes the form: 

(^A, _ ^A,) ^ J2 mv'^^G^aAu, v)-{u^v)]=0. (4.1) 

i 

The sum is over all primary operators Oi appearing in the OPE 0x0, with Aj, Zj, /j 
their dimensions, spins, and OPE coefficients. The unit operator contribution is separated 
explicitly. 

In the approach of [1], one views (4.1) as an equation for infinitely many unknowns 
jf ^ 0. For any given spectrum of Oj's one can ask if a solution exists, and if not, then a 
CFT with such a spectrum is impossible. It turns out that even mild assumptions, like the 
absence of a scalar primary below a certain dimension, can lead to an inconsistent spectrum. 
Making one OPE coefficient too large can also produce an inconsistency. It would take too 
much time to review in detail the results of these studies (see [1, 2, 4-10, 12]), but we would 
like to highlight here two basic issues which were important in all of them. 

First, the functional equation (4.1) looks still rather complicated, and in practice one has 
to replace it by a finite-dimensional constraint, which will of course be weaker but hopefully 
more tractable. One simple way to do it would be to choose a large finite number of points 
in the {u,v) plane and impose (4.1) at each of these points. Ref. [1] followed another way, 
which became standard: to impose (4.1) in Taylor expansion up to a fixed large order 
around just one, well-chosen point. The obvious choice for such a point is 2; = ^ = 1/2, 
i.e. u = V = 1/4, which is invariant under the crossing symmetry transformation z 1 — z. 

The second issue is that the conformal blocks appearing in (4.1) are complicated func- 
tions, even for (1 = 4 where explicit expressions in terms of hypergeometrics are known. Some 
simplifications occur for specific values of A, but this does not help, since the dimensions 
of Oi are unknown and should be allowed to vary freely between the unitarity bound and 
infinity. For generic A, conformal block derivatives at z = z = 1/2 must be evaluated nu- 
merically.^^ For this reason all the studies cited above used numerical analysis. Moreover, 
this evaluation is an expensive operation and often presents a computational bottleneck. 

We will now describe new ways of approaching these issues, made possible by the p-series 
representation. 

4.1 Inexpensive derivative evaluation for all A and / 

Let us briefly review the existing ways of evaluating conformal blocks and their derivatives 
at the point z = z = 1/2. For even d, one uses the explicit representations of Dolan and 

^^For even d, these derivatives can be written via the 3F2 functions (see appendix B.l of [5]). It is not 
known at present how to use these analytic expressions in practice, rather than as a starting point for the 
numerical evaluation. 

^^We should also mention three different types of bootstrap analyses where analytic results could be 
obtained: [3] in the large N expansion; [11] for CFT in presence of a boundary; [13] and [14] in the 
Minkowski space near the light cone. 
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Osborn [19, 20]. For d = 4 they take the form: 



zz 



[/^a+/(^)/;;a-/-2(2;) -{z ^ z)\ . 



(4.2) 



z — z 



By this formula, partial derivatives of Ga.K-^? ^) can be represented as quadratic forms in the 
derivatives of the function ka{z), defined in Eq. (3.22). One can now create an interpolated 
lookup table of ka{z)'s derivatives at z = 1/2 for a range of a. This is a time-consuming 
operation, because the hypergeometric function in (3.22) is expensive to evaluate. However, 
one needs to do this only once. Once the table is created and stored, partial derivatives of 
(4.2) can be computed quickly for any A and /. Such a strategy was used in [2, 4, 6-8], 
and a similar one in [5]. More recently, Ref. [9] found a way to dispense with the lookup 
table altogether, computing the derivatives of ka{z) aX z = 1/2 via a very rapidly convergent 
infinite product representation. 

Turning to general d, a method to evaluate conformal block derivatives was developed 
last year in [10], where it was used to study the 3d Ising model. This method combines 
a variety of ideas. One begins by evaluating partial derivatives along the z = z line, first 
for / = and / = 1 using explicit 3F2 expressions found by [10], then for higher / using 
the recursion relations from [21] reducing those blocks to the lower-spin ones. Then, partial 
derivatives in the orthogonal direction are computed d la Cauchy-Kovalevskaya, using the 
fact that the conformal blocks satisfy a second-order partial differential equation. 

The p-series gives a new way to evaluate conformal blocks and their derivatives, which 
works for general d and around any z. To achieve the necessary precision, one needs to 
evaluate the coefficients Bnj as a function of A and Z up to a sufficiently high order, using 
the recursion relation (3.12). It is important that the necessary number of terms will be 
independent of A and I, because of the bound (3.18). For example, to be able to compute 
the conformal blocks with double precision (10^^^) one would need the coefficients up to 
level 



which gives n ^ 20 for z = 1/2. This number is a bit of an underestimate, because it 
assumes that the sum of the coefficients at level n is uniformly bounded, while in fact it 
grows with n as in Eq. (3.21). Also more levels will be needed if one wants to evaluate 
derivatives. 

It should also be rather easy to generalize the p-series method to the case of unequal 
external dimensions Ai 7^ A2 and A3 7^ A4. The extra terms in the Casimir operator for 
unequal external dimensions are all first order in derivatives [20]. So, the operator V will 
remain first order, and we can expect that the boundedness properties of the coefficients 
Bn,j will still hold. Such a generalization will be useful for the conformal bootstrap analysis 
of several scalar correlators simultaneously. For example, in the 3d Ising model, it would be 
interesting to study simultaneously the correlators {aaaa), {aaee) and (eeee). The second of 
these correlators has unequal external dimensions in two conformal partial wave expansion 
channels out of three. 



16/logio(l/p), 



(4.3) 
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4.2 Analytic toy model for the conformal bootstrap 

As already mentioned, Eq. (4.1) is usually analyzed numerically. Numerics are essential both 
for computing the conformal blocks and their derivatives at z = 1/2, and for performing 
searches in the resulting derivative spaces, which in the most advanced studies [9, 10] can 
have O(IOO) dimensions. If one needs precision, numerics are unavoidable at present. 

Leaving precision aside, here we would like to address a more modest question: can one 
provide an analytic understanding of why the method of [1] gives nontrivial constraints? 
Some intuitive explanations were given in section 5.1 of [1], but those still relied on properties 
of conformal blocks which had to be checked by plotting them in Mathematica. As we will 
now explain, analytic results can be obtained using the p-series representation. 

The basic idea is the following. We have shown that the conformal blocks can be repre- 
sented as series in the "Gegenbauer blocks" ^^^(r, 77). These series are rapidly convergent, 
to the extent that even the first term provides already a very good approximation: 

GA,i^VA,i{r,v)- (4.4) 

The relative error is of order p^, which is about 3% aX z = 1/2. Importantly, the error 
is uniformly small for all A and /. Replacing the conformal blocks in Eq. (4.1) by their 
Gegenbauer block approximation, we get the "toy" bootstrap equation {5 = A^): 

[(1 - zY{l - zY - /z'] + ^ IIiIHaAz, z) - i^A,/(l - ^, 1 - ^)] = , 

H^A^rz) - (1 - - -zY[{p{z)p{-zri'c\^' ( 2[p(1)p(S/0 • ^^-^^ 

It is expected to give qualitatively the same results as the full bootstrap equation, with an 
advantage that the analysis can be done analytically, since the Gegenbauer block derivatives 
can be computed explicitly. We will give below two examples of this approach. 

A comment is in order concerning the spectrum of A's and /'s appearing in the second 
term of the toy bootstrap equation. If we view it as an approximation to the full equation 
(4.1), then of course it's the same spectrum as in the full equation, i.e. all primaries in the 
0x0 OPE. However, an alternative point of view can be useful. We can consider the toy 
equation as exact, provided that we allow not only the primaries but also their descendants 
to appear in the spectrum. When we do the bootstrap in terms of the full conformal 
blocks, we have extra constraining power because the OPE coefficients of descendants are 
proportional to those of the primaries. In the toy bootstrap, we choose to discard this 
information and allow descendants to appear with independent coefficients. One can also 
imagine an intermediate situation, when descendants up to a certain level are included 
with the relative coefficients fixed by conformal symmetry, while the higher ones are taken 
independent. 
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4.2.1 Toy bootstrap for z = z 

For the first example [32], let us consider the toy equation (4.5) for real < 2; < 1. The 
angular part of the Gegenbauer blocks is then trivial, and the equation takes an extremely 
simple form: 

[(1 - zf - z''] + /!{(! - zr[p{z)]'^ - z''[p{l - z)]^} = . (4.6) 

A 

We will use this equation to show that there is an upper bound on the lowest primary 
dimension Amin in the 0x0 OPE. This is a problem of the kind first considered in [1], 
except that here we are not distinguishing between scalar and higher spin primaries. 

For the proof, let us Taylor expand Eq. (4.6) m x = z — 1/2. Only odd powers of x will 
appear since the functions are odd. From the first term we get: 

(1 _ _ ^ (x + 1(5 - 1)(25 - l)x3 + ...), (4.7) 

where > is a constant whose precise value is unimportant. When expanding the 
Gegenbauer block terms, let us assume that all A ^ 5 (we aim for a contradiction here). 
Then we can approximate z'^^ ~ (1 — zY^ ^ (1/2)^*^ = const, while the relevant part is: 

[p{z)f - [p(l - z)f = B^{x + lA'x' + ...) , (4.8) 

where B/\ > is another inessential constant. We can change normalization of the blocks so 
that Ba 1 (incorporating this constant into /^). Requiring that (4.6) be satisfied term 
by term in the Taylor expansion, we get: 

= E /a , CsiS - 1)(2<5 -l) = Yl ^'fl ■ 

The 0{x^) terms etc. would give more equations but we won't use them here. Bounding 
the RHS of the second equation from below by A^j^ ^ and using the first equation, we 
conclude that 

A^in ^ V(5-l)(25-l) . (4.9) 

This shows that our original assumption that Amin ^ 5 is inconsistent, hence there must 
be a bound on Amin in terms of S. Its actual value can be found by a more careful analysis, 
expanding the Gegenbauer block terms without the approximation A ^ 6. 

4.2.2 Including the spin dependence 

We will next consider the toy bootstrap equation (4.5) not restricting to the z = z line. 
The advantage is that the spin information will be now accessible through the order of the 
Gegenbauer polynomial. So we can try to set a upper bound on the lowest scalar in the 
(j) X (f) OPE, which is precisely the problem considered in [1]. 

For the simplest bound, we will expand (4.5) to the third order in z and z around 
z = z = 1/2. Because of various (anti) symmetries, only three derivatives are independent. 
We will choose 

9„ dl-d^d,, dl + 3d% (4.10) 
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as a basis, as these linear combinations somewhat simphfy the subsequent algebra. Because 
all functions are elementary, the derivatives can be evaluated explicitly. Equating them to 
zero, we get the following linear system: 

5^p(l/2)^/i,hA,, = ho,o, (4.11) 

A,l 

ho,o = [25, -16(5(5 - 1), 165(2(5^ -35 + 1)]* , 

hA,i = [v^A - 25, -4A2 + A V2(-45 + Sk + 7) + 165^ - 165{ki + 1) - 8ki, 
8V2A^ - 12(45 + 1)A2 + ^2(4852 - 125 + 7)A - 165(25^ -35 + 1)]* , 

where ki is the logarithmic derivative of the Gegenbauer polynomials at r/ = 1: 

lc!+mjl±2^^ (4,12) 

Notice that the spin and the spacetime dimension enter only through this coefficient and 
via the unitarity bounds. 

Next we eliminate the RHS from the last two equations in (4.11) with the help of the 
first one. We get a homogeneous system: 

J2 ^A,/ gA,; = (0, 0)* , qA,i = AV(l/2)^/l,, , (4.13) 

1 - 3V2{5 + l/4)A-i + (452 + 35/2 - l/8)A-2 
A-i - V2{5 + 2ki - l/4)A-2 + (45 + 2)A;;A-3 



(4.14) 



The rest of the discussion follows closely section 5.4 of [1]. We have to study how the 
direction of the vectors gA,; varies when we increase A from the unitarity bound to infinity. 
For / = we vary A from Ao,min to infinity, where Ao,min is the lowest scalar dimensions. If 
the set of all directions stays within a cone of opening angle < tt, then (4.13) will not have 
a nontrivial solution, while in the opposite case it will. The difference from [1] is that now 
Sa,i are given explicitly, and the analysis can be carried out analytically. We will not give 
here the details, but we have checked that an upper bound on Ao,min can be obtained using 
this method for all 2 ^ (i ^ 4, at least for 5 near the scalar unitarity bound. 



4.3 Truncated bootstrap equation with an error estimate 

As mentioned above, one could also try to do bootstrap imposing the bootstrap equation 
point by point at several z = Zi, rather than in the Taylor expansion around z = 1/2. 
We would like to discuss here the issues arising if one wants to implement this technique. 
Conformal block evaluation for any z can be done with the p-series. The next question 
is then how to distribute the sampling points. To get an idea, let us consider the rate of 
convergence of the conformal block decomposition (2.3). As shown in Ref. [27], the error 
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induced by truncating (2.3) at some maximal dimension A = A,,, is exponentially small: 



< 



r(4A^ + i) 



(4.15) 



To be precise, this estimate was shown to hold for A^, ^ A0/(1 — |p(-2) |) . Most importantly, 
it holds in an arbitrary CFT with no extra assumptions about the 0x0 OPE. For example, 
it might seem that having too many operators at high A, or a single operator with a huge 
OPE coefficient, might invalidate this bound. However, the proof in [27] shows that such 
situations cannot occur in a consistent CFT. 

The estimate (4.15) is relevant to our discussion, because in most practical approaches to 
the bootstrap one has to truncate the spectrum of considered operators from above (to make 
the problem finite). Now we know that the error induced by this operation is controlled by 
|p(^;)|, while the error in the crossed channel will be controlled by |p(l — z)\. Therefore it 
seems natural to distribute the points Zi in a region of the form (see Fig. 9) 



X{z) = max(|p(2;)|, |p(l - z)|) ^ , 
where Ac should be chosen commensurately with the eventual dimension cutoff A^,. 

2.0- 



(4.16) 




Figure 9: The contour plot of the function X{z) in the plane (Re z,lm z). Only the region 
Im z ^ is shown, since the conformal blocks are symmetric in z, z. 



One way to choose A* is so that the error (4.15) is below the numerical precision one is 
working with (say double precision) everywhere within the region (4.16). Alternatively, one 
can choose A^, lower, so that the error is non negligible. Then one has to include this error 
estimate directly into the bootstrap equation. Such a modified equation takes the form: 



u 



A(C')s£A. 



^ Eiz^z) 



4Ai 



Eiz.z) 



A 

r(4A^ + 1) 



max 



|1 - z\^^^\p{z)\^% kl'^1p(l - ^)|^*) • (4.17) 
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We think it would be interesting to try to carry out bootstrap analysis based on this 
"truncated bootstrap equation" rather on the conventional technique of expanding around 
z = 1/2. There are many free parameters one can play with: A*, Ac, the number of points 
Zi at which to impose (4.17), and what is the optimal way to distribute them in the region 
(4.16). Once all these parameters are fixed, the problem of deciding whether (4.17) has a 
solution with ^ can be solved via the linear programming algorithms. 

It's worth pointing out an additional feature of Eq. (4.17), which makes it particularly 
useful when the conformal blocks are computed via the p-series, whose coefficients can be 
computed up to arbitrary order but whose closed form is unknown. Namely, it remains valid 
when the conformal blocks Go{u,v) are replaced by the "truncated blocks" — the partial 
sums of the p-series up to the level A + n ^ A^. This is because the error estimate (4.15) 
is in fact valid when the contributions of all states of dimension above A.,, are included into 
the LHS (and not just the conformal multiplets of primaries above A.,,). It is in this stronger 
form that the error estimate was proved in Ref. [27]. 

5 Discussion 

In this paper we developed the theory of conformal blocks rooted in their physical meaning 
as sums of exchanges of descendant states in the radial quantization. This point of view is 
standard in the 2d CFT literature [26], and our goal here was to demonstrate its utility in 
higher dimensions. 

We explained how quantum mechanics fixes the structure of conformal block in radial 
coordinates: it is an integer-spaced power series in r with angular dependence given by 
Gegenbauer polynomials. The coefficient of each term is positive as a consequence of 
unitarity. These coefficients are easy to find using recursion relations following from the 
fact that the conformal blocks are eigenfunctions of the quadratic Casimir of the conformal 
group. 

We highlighted the existing freedom in the choice of the radial coordinates. It's the same 
freedom as when expanding the product of two operators 0(xi)0(x2) into a sum of operators 
inserted in some point Xq, which becomes the radial quantization origin. Each choice gives 
a different representation of the same conformal block, and it is not a priori clear which one 
is more convenient. In this paper we analyzed in detail two natural choices, the end point 
xq = xi, and the middle point xq = (xi + X2)/2. 

The end point choice (section 2) corresponds to working with the complex variable z 
often used to represent conformal blocks, with explicit 2-^1 representations available in even 
dimensions d. For general d considered here, we expand conformal blocks in a power series 
in \z\ times Gegenbauers. The expansion coefficients Anj satisfy a three-term recursion 
relation, also derived earlier from a different point of view by Dolan and Osborn [20]. An 
unpleasant feature of these expansions is that the coefficients at level n grow with the 
exchanged primary dimension as A". For large A many terms need to be evaluated to get 
a good approximation to the conformal block. 
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Choosing the middle point (section 3), one passes from z to the complex variable 



P = . r. ■ (5-1) 



z 



This variable was recently used in Ref. [27] to study convergence of the conformal block 
decomposition. As we showed here, this is also an ideal variable for constructing rapidly 
convergent expansions of the conformal blocks themselves. The expansion coefficients 
satisfy a recursion relation which is a bit more involved than for the v4„ the coefficients 
at level n are linear combinations of coefficients at all levels n — 2, n — 4, ... up to zero, 
while for An^j only the level n — 1 contributes. But this complication pays off: the resulting 
coefficients do not exhibit any growth with A or /. This means that the coefficients computed 
and stored up to some large and fixed level N can be used to evaluate conformal blocks of 
arbitrary dimension and spin with uniform accuracy. 

Amazingly, even the first term in these infinite p-series expansions provides already a 
pretty good approximation (within a few %) to the full conformal block. We use this fact in 
section 4.2 to propose the "toy bootstrap equation". Although this equation discards some 
information compared to using the full conformal blocks, and hence is less constraining, 
it has an advantage of involving only elementary function and being amenable to analytic 
analysis. Using this "toy bootstrap" , we get analytic understanding of why the method of 
[1] was able to get an upper bound on the lowest dimension in the OPE. 

We believe that our p-series expansions will find many other future uses in the bootstrap 
program. Some of the possibilities are described in section 4. We conclude here by 
considering another possible application: bootstrap analysis of four point functions of non- 
scalar external primaries. Conformal blocks for such correlators have been studied recently 
in [22, 23]^°, but the results were not yet put to concrete use. Partly this is due to the 
fact that the obtained expressions are still rather complicated and not fully general. For 
example, Ref. [22] finds only the blocks corresponding to the symmetric traceless exchanged 
primaries, while more general representations can be exchanged if the external fields have 
spin. We believe that the p-series approach could be useful in the problem of expressing these 
missing conformal blocks. The basic building blocks will no longer be simple Gegenbauers, 
but they will still be polynomials of the angular variable, fixed by the SO{d) group theory. 
Once the expansion basis is known, the coefficients can presumably be found by using the 
Casimir equation judiciously. It would be interesting to carry out this computation in detail. 
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A Boundedness of the p-series coefficients 



In this appendix we show that the coefficients Bnj on each level n are uniformly bounded 
for all A and I, as stated in Eq. (3.18). We have already shown in the main text that Bnj 
remain bounded as A — )■ cxd for each fixed I. So here it suffices to consider the case of I large 
with respect to n, say I ^ n. 

We will proceed by induction, and assume that the inequality has already been shown 
for all levels n' < n. Using Eq. (3.12), the bound at level n will follow if we show that 

^iTn',' /(CA^n, - C^,l) (A.l) 

is bounded by a constant which depends only on n and u. 
Our ffist observation is that the F's satisfy the bound 

^A+n'''j' ^ const. A + const (A. 2) 

with constants which depends only on n and u. To show this, notice that large contributions 
to F's appear from only two sources. First, through the action of rdr, which gives a factor 
of (A + n'). Second, through the action of 27^(1 — 77^)^^, which gives factors bji = 0{j'). 
On the other hand, all the factors produced via the expansion of denominators in (3.7) will 
depend only on n and u. Notice in particular that aj = 0(1). 

Passing to the Casimir difference in (A.l), we write it as 

Ca+uj - Ca,i = k'^ + 2k{l - n + z/) + 2n(r + n - 1) , (A. 3) 

where j = I — n + k, k = 0,2, ... , 2n, and r = A — / — 2z/^0by the unitarity bounds. Since 
we are assuming / ^ n, this is a manifestly monotonically increasing function of k and n. 

Consider ffist the case k ^ 2. In this case we have a lower bound: 

CA+nj - Ca,i ^ [(A.3) for = n = 2] = 4(A - u) {k ^ 2) . (A.4) 

Combining this with (A. 2), we see that (A.l) is indeed bounded independently of A and I, 
except in the region near the free scalar unitarity bound A = u, excluded from consideration 
as discussed in the main text. 

It remains to consider the case k = 0, when the Casimir difference 

CA+n,i-n - Ca,i = 2n(r + n - 1) (A.5) 

can remain small even though both A and / become large. However, precisely in this case 
the bound (A. 2) can also be improved. The relevant recursion relation coefficients are 

^"^J^i-w ■ (A.6) 

The coefficients with j' I — n' will be zero, because lowering the spin via Eqs. (3.9) is 
accompanied by raising the dimension by at least the same amount. We will now show that 
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all coefficients of the form (A. 6) are bounded by const. t + const. Together with (A. 5), this 
will prove the boundedness of Bn,i-n and will complete the proof. 

For the case n — n' = 2, this stronger bound can already be suspected in the expression 
(3.11) for r^^^'"'"^, which contains two near-canceling terms. In detail, this coefficient can 
be expressed as 

rS::5!l7"'-^ = 4(r + 2n')ar.„, , (A.7) 
and satisfies the claimed bound, since aj ^ 1/2 for all j. 

For the general case, we notice that the action of V on VA+n',i-n' can be written as 
follows: 

A{t + 2n')aY_^, 

VVA+n',l-n' = - ^ _ 2^2(27^2 _ ^ ) ^A+n'+2,«-n'-2 + • • • (A.8) 

Here we computed explicitly the action of (1 — 2r]'^)rdr and 2ri{l — i]'^)drj. We omitted many 
terms (...) which cannot contribute to the relevant F coefficients, because they raise the 
dimension without lowering the spin by the same amount. The F coefficients (A. 6) with 
n — n' = 4,6, . . . are obtained by expanding the denominator in (A.8). They are given by 

riStn' = 2(r + 2n') J] (^aj) , (A.9) 

j=l-n',l-n'-2,...,l-n+2 

and clearly satisfy the claimed bound. 

The reader will have noticed that the coefficients Bn^i-n, which required a separate 
analysis in the above proof, satisfy a recursion relation among themselves. This is because 
the relevant F's in (A. 6) vanish for j' ^ I — n' . Due to this fact, these coefficients can in 
fact be computed exphcitly: 

_ (l/2)^ {l + l-2m)2m {r/2)m 
^2m,i-2m ^, ^ ^ + ^ _ 2m)2^ (r/2 + 1/2)^ ' ^'"•''^^ 

Their boundedness is also obvious from this formula. 
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